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We microscopically model the decoherence dynamics of entangled coherent states of two 
optical modes under the influence of vacuum fluctuation. We derive an exact master 
equation with time-dependent coefficients reflecting the memory effect of the environ- 
ment, by using the Feynman- Vernon influence functional theory in the coherent-state 
representation. Under the Markovian approximation, our master equation recovers 
the widely used Lindblad equation in quantum optics. We then investigate the non- 
Markovian entanglement dynamics of the two-mode entangled coherent states under 
vacuum fluctuation. Compared with the results in Markovian limit, it shows that the 
non-Markovian effect enhances the disentanglement to the initially entangled coherent 
state. Our analysis also shows that the decoherence behaviors of the entangled coherent 
states depend on the symmetrical properties of the entangled coherent states as well as 
the couplings between the optical fields and the environment. 
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1 Introduction 

Optical fields are widely used in quantum communication since quantum information is almost 
invariably transmitted using photons. Experimental quantum teleportation has been realized 
using the discrete two-photon polarization entanglement states or the continuous two- 
mode squeezed entanglement states [2j |3] as quantum channels [J . Another important type 
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of continuous variable entanglement states, entangled coherent states O El [7l [3 [9j [10], has 
also been proposed as a potential quantum channel to teleport unknown quantum states 
[moils]. In this paper, we shall investigate the non-Markovian decoherence dynamics of 
the continuous variable quantum channel in terms of entangled coherent states. 

As well known, a realistic analysis of quantum systems for quantum information pro- 
cessing must take into account decoherence effect. There has been an increasing interest in 
describing various continuous variable quantum channel under noise [Ml [TS] [TBI [IZl [Ml [H] • 
Conventional approaches treat the interaction between the quantum system and its environ- 
ment perturbatively, which yield equations of motion such as Redfield or master equations 
under the Born-Markov approximation [2Ql [211 [22] • Although the approximation has been 
widely employed in the field of quantum optics, where the characteristic time of the envi- 
ronmental correlation function is short compared with that of the system investigated [22j . 
its validity is experiencing more and more challenges in facing new experimental evidences 
[23j . Moreover, the Born-Markov approximation is in general invalid in dealing with most 
condensed-matter problems, for example, a quantum system hosted in a nanostructured en- 
vironment [24l |25l [28l |26l [27] , because possible large coupling constants and long correlation 
time scales of the environment both require a non-perturbative treatment. Thus, how to 
develop a general non-perturbative microscopic description of open quantum systems has 
attracted much attention recently [29l [30l [3ll [32] . 

In the present work, we shall focus attention on the influence of vacuum fluctuation on 
quantum channels in terms of entangled coherent states. To this end, we model the system 
as two optical modes coupled to a bosonic environment at zero temperature. We shall then 
develop non-perturbatively a microscopic description to the decoherence dynamics of such 
systems. We have noticed that most of previous theoretical works to explore the decoherence 
dynamics of the optical field system relied on Born and/or Markov approximation [Ml [151 
[T6l \YJ\ [TSl [T9] . To derive non-perturbatively the decoherence dynamics of an open quantum 
system, we will employ the Feynman- Vernon influence functional theory [33l [HI [35] in the 
coherent state path integral formalism [36j , which enables us to treat both of the back-actions 
from the environment to the system and from the system to the environment self-consistently. 
After a careful evaluation of the coherent state path integrals, we obtain an operator form of 
the exact master equation with time-dependent coefficients describing the full non-Markovian 
dynamics of the back-actions between the system and the environment. 

We then investigate the non-Markovian decoherence dynamics of the entangled coherent 
states [5] using the exact solution of the reduced density matrix, where the entanglement 
is measured by the concurrence [37]. The non-Markovian effect is manifested in the short 
time peak of the time dependent coefficients in the master equation, which results in an 
enhancement of the disentanglement to the entangled coherent states. Indeed, in a recently 
published paper |32j , we have already used the exact non-Markovian master equation derived 
in this paper to study the decoherence dynamics of another type of the continuous variable 
entangled states, i.e. entangled squeezed states, where the entanglement is determined by the 
logarithmic negativity |38j rather than the concurrence since the later is not applicable to 
entangled squeezed states. 

The paper is organized as follows. In Sec. II, we introduce the model of two optical modes 
interacting with a common environment in the coherent-state representation. In Sec. Ill, we 
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show the detailed derivation of the influence functional theory to the model. The exact master 
equation is derived in Sec. IV. Sec. V is devoted to the study of entanglement dynamics and 
the decoherent properties of quantum channels in terms of the entangled coherent states. 
Finally, a brief summary is made in Sec. VI. 

2 The Hamiltonian of two optical modes in an environment 

Our system includes two separated optical modes subject to a common vacuum fluctuation, 
which is relevant to quantum network and has been widely investigated j391 140j . Since we are 
interested in the decoherence of the two optical modes mediated by a vacuum electromagnetic 
field after the two-mode entangled coherent state is prepared [71 IE] , we can omit the terms 
regarding the atoms in 30] . The Hamiltonian of the whole system is then given by 

H ^Hs + He + Hi, (1) 

where 

Hs — hhJia\ai + fiu}2a\a2 + 'hK{a\a2 + a\ai), (2) 

He = ^^TiuJkb\bk, (3) 

k 

Hi = J2 Hgikajbk + gtk^ibl), (4) 

l,k 

are, respectively, the Hamiltonians of the two optical modes, the environment (vacuum fluc- 
tuation), and the interaction between them. The operators a/ and a| {I = 1,2) are the 
corresponding annihilation and creation operators of the Z-th optical mode with frequency lui . 
The parameter k is a coherent tunneling rate of photons between the two optical systems, 
such as two cavities [4TJ |42] , which is proportional to the overlap of the two wave packets 
of the optical fields. Such coupled optical array system recently attracts much attention 
[ITl [53] for the possible materialized in a variety of physical systems, for example, fiber 
coupled micro-toroidal cavities [H] , arrays of defects in photonic band gap materials [IS] and 
superconducting qubits coupled through microwave stripline resonators |46j . 

The environment is modeled, as usual, by a set of harmonic oscillators identifying the 
vacuum electromagnetic field with the annihilation and creation operators bk and bl.{k = 
1, 2, • • • ), gik are the coupling constants between the optical modes and the environment. In 
Eq. ([T]) we have also suppressed the polarization of the fields for both the systems and the 
environment. Since most quantum optical experiments are made currently in low temperature 
and under vacuum condition, the vacuum fluctuation should be a main source of decoherence. 
Therefore, we take the environment to be at zero temperature throughout this paper. 

To apply the influence functional method to an open quantum system, the first step 
towards the dynamics of the reduced system is to compute the forward and backward prop- 
agators between certain initial and final states of the full system by choosing a convenient 
representation. In the present work we use the coherent-state representation [36], in which 
the basis of the Hilbert space for the environment consists of multi-mode bosonic coherent 
states 

|z) = n l^*^)' 1^'=) " exp{zkbl)\Ok), (5) 

k 
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and that for the two optical modes is the two single-mode bosonic coherent states 
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|Q!) = ni"'''' I"') = exp(a/a3')|0,), (6) 



1=1 

where the shortened notations for the complex variables, z = (zi, Z2, • ■ ■ ) and ol = (ai, a2) , 
are introduced. 

The coherent states defined above are eigenstates of annihilation operators, 

bk\zk) ^ Zk\zk) , ai\ai) ^ ai\ai) . (7) 
As these coherent states are over-complete, they obey the resolution of identity, 

j d/.(z)|z)(z| = 1, y" d^l{a)\a.){oc\ = 1, (8) 

where the integration measures are defined by c?/i(z) = e~^'=^'= '^^2*th''° ^'^'^ '^f^ ('^) ~ 11/ '~2^P'- 
As it is shown, the bosonic coherent states we used here are not normalized, and the nor- 
malization factors are moved into the above integration measures, which corresponds to the 
Bargmann representation of the complex space. Moreover, these coherent states are also 
nonorthogonal, 

(z|z') = cxp(^ zlz'k), {a\a') = exp(^ a^a'i). (9) 

k I 

The use of the coherent-state representation makes the evaluation of path integrals extremely 
simple. In the coherent-state representation, the Hamiltonians of the two optical modes, the 
environment (vacuum fluctuation), and the interaction between them are expressed as 



Hs{ot, a) — h ujiaiai + hK{aia2 + a2ai), (10) 
1=1 

He{z,z) =^iUkZkZk, (11) 
fe 

Hi{a, a,z,z) =^{gikaizk + gtkZkai), (12) 

Ik 

where z and a denote the complex conjugate of z and a, respectively. With the above 
coherent-state representation, we will present in the next two sections a detailed derivation 
of the exact master equation for the reduced density matrix of the two optical fields that we 
have simply outlined in our early work |32j . 

3 The influence functional theory 

3.1 The influence functional in coherent- state representation 

We follow the influence functional method of 47J by expressing the density matrix of the 
composite system as a double-path coherent state path integral. After eliminating the degrees 
of freedom of the environment, we can incorporate all the environmental effects on the reduced 
system in a functional integral named influence functional [33] . Then the dynamics of the 
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reduced system will be governed by an effective action retaining all the influences from the 
environment in the influence functional. 

The total density matrix of the system plus the environment obeys the quantum mechan- 
ical equation ihdptot(t)/dt = [77, ptot(i)], which yields the formal solution: 

Ptotit) = e^ptotiO)e'^. (13) 

Different from the coordinate representation in [34l [35], the coherent-state representation 
leads to, 

(a/,z/|ptot (t) \a'f,Zf) 

= J dp{zi)dp{a.i)dp{z[)dp{a'.i){af,Zf:t\oti,Zi]0) 

X ( Q;j,z,|ptot(0)|a-,z',)(Q:^,z^;0|a'y,Z/;i}, (14) 

where the resolutions of identity, Eq.®, has been used. The density matrix given by Eq. 
describes the behavior of the two optical modes plus the environment as a whole. As we 
are only interested in dynamics of the two optical modes, we will work with the reduced 
density matrix by integrating over the environmental variables. We also assume that the 
initial density matrix could be factorized into a direct product of the two-mode state and 
the environment state /f)tot(O) = p{0) (8> P_e(0), namely, assuming no correlation between the 
environment and the system at t < |48j . Then the reduced density matrix fully describing 
the dynamics of the two optical modes is given by 

p{af, a'j; t) ^ J dp,{oLi)dp{a'i)p{a.i, a^; 0) 

X J{a.f,a.'f,t\a.i,a.[]Q), (15) 
where p(q:,q:';t) = J d/x(z)(Q:, z|/9tot(''")|cK', z), and 
J{a.} ,a'f\t\a.i,a[\Q) 

= j dp{zf)dp,i7.i)dp.{z[){cy.f,Zf:t\a.i,Zi\Q) 

X P£;(z»,Zi;0)(Q!-,z^;0|aj,z/;i), (16) 

is the propagating function of the reduced density matrix, which contains two propagators 
for the total system; the forward and backward propagators, e^^, plus the initial density 
matrix of the environment as a matrix element in the coherent-state representation. 

In the following we will show how to calculate the forward propagator in terms of the 
coherent state path integral [49l [36] . The similar calculation could be done for the backward 
one. To evaluate the forward propagator operator e k between the initial {\a.i^Zi)) and 
the final ((a/,zj|) coherent states, one can generally divide the time interval tf — ti into 
N subintervals. Then by inserting the resolution of identity {N — 1) times between each 
subintervals and taking the limit of large N , we have the forward propagator in terms of the 
coherent state path integral, 

{a.f,Z}]t\ ai,Zj;0)= j ZJ^zZ)^ a exp |^ (55[q:, a] 

+ 5/[z,z,a,Q;] + S'£;[z,z])}, (17) 
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where Ss, Se, and Sj are the actions corresponding to the two optical modes, the environment, 
and the interaction Hamiltonian iJg, He, and Hi, respectively, 

Ss[ot,ot] = X] { ^ T-^^io^i (t) 

+ J dT[ihaiai{T) ~ Hs{a,a.)]Y 
Se[z,z] =X] { ~ iTT-ZkZkit) 

k 

+ j dT[i%ZkZk{T) - He{z,z)\^, 

^/[z, z, Q:,a] = — / dTHj{a., a.,z,z). 
Jo 



(18) 

(19) 
(20) 



All the functional integrations are carried out over paths z(r), z(t), q:(t), and q:(t) with 
endpoints z{t) = zj, z(0) = Zj, a{t) = a/, and a(0) = a^. Substituting Eq. (fT7|) and a 
similar expression for the backward propagator into Eq. (|16p , we obtain 

J'{af,a'j;t\ai,a^;0) = J D^aD^a' cxp ^^{Ss[a, a] 

~Ss[a', a'])^T[a, a, a', a'], (21) 

where 



T[a, Oi, a' , a'] = J dii{zf )dii{zi)dii{z^)D^zD^z' 

X P£;(Zi,Zi;0)exp|^(S'B[z,z] ~ 5|;[z',z'] 

+ S'/[z,z,a,a] - 5;[z',z',a',a'])| (22) 

is defined as the Feynman- Vernon influence functional in the coherent state representation, 
which contains all the environmental effects on the two optical modes. 

3.2 Evaluation of the influence functional 

Now we can calculate explicitly the influence functional of our model using the coherent-state 
path-integral formalism presented above. Substituting the model Hamiltonian into the actions 
of Eq. p8ll20p . we obtain the explicit form of the forward propagator. The path integral of 
the environmental part of the propagator can be done by the stationary phase method [491 136] 
with the boundary conditions 2:^(0) = Zki and Zk{t) = Zkf, which results in the equations of 
motion, 

Zk + ioJkZk = -i^^gtkCti, Zk - iuJkZk = i'^gikOii, (23) 
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where the paths regarding ol and a are taken as external sources. The solution to the 
stationary path equation ([23]) are 



, J T 



(25) 



Note that the prefactor under the contribution of stationary path in the coherent-state path 
integral is unity, and the stationary phase method to treat the environmental part here is 
exact for the action being only a quadratic function of the dynamical variables. The path 
integral of the environmental part for the backward propagator (a^, z^; Ojap zy; can be 
evaluated in the same way. 

Since we only consider the vacuum fluctuation, the environment is initially in the equi- 
librium state at zero temperature, we then have /5£;(zi, z^; 0) = 1. Substituting the solution 
(I24II25P for ^^(t), Zfc(T) and a similar solution for z'f.{T), z'i^{t) together into Eq. (|22p . and 
using the Gaussian integral identity / d_^^~izz+\z+}yz _ ig— repeatedly for the integral over 
Zi,z[,Zf, we reach the final form of the influence functional that we have used in [32], 

J-[(y, a, a' , a'] = 

expj / dr dr' V(a;(r) - ai(T))/Xi„(r - r') 
X a™(r') + {ai{T) - a[{T)) ^l*„^{T - /^(t')] }, (26) 

where fiim{x) — '}lk^^^'^''^ dikdmk dissipation-noise kernel. 

4 The exact non-Markovian master equation 

4-1 The propagating function of the reduced density matrix 

In the above derivation of the influence functional, the back-actions between the two optical 
modes and the environment have been treated self-consistently. All the effects from environ- 
ment on the two optical modes are incorporated in the influence functional which leads to a 
modification to the action of the two optical modes, 

J{a.f.,OLj ]t\a.i.,a.[]0)= I D^olD^ol exp | {aiai (t) 

/■* ^ _ . 
+a'ia'i{t)) - / dT[y2{aiai + a'la'i) + iHs{a,a) 

•^0 1=1 

—iHs{a.',a')^^!F[6i, a, a', a']. (27) 
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To execute the path integral of Eq. (|27p , again we resort to the stationary phase method and 
obtain the equations of motion as [l ^ V ) 

ai + iiuJiai + kq;//) = - / dr' ^ (t - r') (t') , (28) 

2 

a[ - z{LOia[ + K«;,) = - f dr' ^ (t - r') a'^ (r') . (29) 

with the boundary conditions ai (0) = an and (0) = aj^. 

The integro-differential equations render the reduced dynamics non-Markovian, with the 
memory of the environment's dynamics registered in the time-nonlocal kernels. To simplify 
the discussion, we further assume that the two optical modes are identical, i.e., wi = = cjq. 
Then the coupling strength to the common environment should also be the same: gik ~ 
e"^(72fc = 9k, where the phase factor e'*^ = A models the phase difference between the two 
optical modes coupling with the environment. In the present work we will consider two 
special cases: the two optical modes couple with the environment in phase (a constructive 
interference coupling with ip — X = 1) and out of phase (a destructive interference 
coupling with = tt so that A = —1). By introducing the new variables 

ai{T) = auuir) - ai>iv{T), 

a[{T)=a[,u{T)-a[,^v{T), (30) 

and using the equations of motion (|28ll29p , we obtain the propagating function of the reduced 
density matrix as 

J{a.f,a'j;t\6n,a'f,Q) = 
2 

exp I \uaifaii + ua'na'i^ — {uu + vv — l)ciij^aij\ 
1=1 
2 

- ^ [vaifai>i + va'^a'i, ^ - {uv + vu)a'iiai'i] |. (31) 
{1,1') 

where u, v are solutions of Eq. pop at time t = t. The exact reduced density matrix is 
then easy to be obtained by substituting the above solution of the propagating function into 
Eq. (fT5)) and integrating over the initial state. 

4-2 The exact non-Markovian master equation 

Eq. ([3T|l is an exact result. In this section, we will deduce the master equation from Eqs. (fT5|) 
and ([31]). From Eq.(l3l|), we obtain 

auJ - 5 5 , Ui^J ~ , (62) 



which will be used to eliminate the dependence on the initial values ai,a^ in Eq. (fT5|) . 
Combining Eqs. ((3T|) and ([15]) together, and using the identities of Eq. ((32| . the evolution 
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equation of the reduced density matrix is given by 



p(a, a ; i) - }^ I - in{t) [a, — — ai\ 

1=1 

SaiSai Sai Sai ' 



{IM) 



dp{a,a';t) 6p{a,a';t) 
ai J ai> 



, T-'^+^ro'^^'°("'"''^) - 5p{a,a'-t) 5p{a.,a.':t) 
oaioai oai oai 



(33) 



where 



T,/ X ^/ N uu — vv 

r{t) + in{t)^ — 2 — 

. N „// X vii — uv , ^ 

T'it)+tn'it) = -- J. 34 

Eq. p3p is the exact master equation of the reduced density matrix for the dynamics of the two 
optical modes in the coherent-state representation, in which fl{t) plays the role of a shifted 
time-dependent frequency of the two modes, Q'(t) accounts for a shifted time-dependent 
coherent interaction between the two modes, T{t) represents a time-dependent individual 
decay rate of each mode, and T'{t) is for a correlated decay rate of the two modes induced by 
the environment. 

If we define a new variable F±{t) = u{t) ± v{t), then Eqs. p8ll29p is reduced to 

+ (1 T A) r dT'p{T - r')F±(r') = 0, (35) 
Jo 

with p{x) = Efce-*"'="^|5fc|2 and A = ±1. The explicit forms of n{t), n'{t), and r(t) in the 
master equation are given by 

il{t) =a;o + Ini[GA(t)], n' (t) = k + Mm [Gx{t)] , 

r{t)^Xr'{t)^Rc[Gx{t)], (36) 



where 



1 



F-xit) 

dTp{t-T)F^x{T). (37) 



This result has the similar form as the coefficients in the non-Markovian master equation of 
a two-level atom derived in [29 . 
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To obtain the operator form of the master equation, we should introduce the foUowing 
functional differential relations in the coherent-state representation (i.e., the Bargniann rep- 
resentation of operators [5Tj). 

- 5p{a., a'-t) I 

Oil FT < mlarnPyt): 

bam. 



5p{cx^ a ; t) 

J dm 

Sai 
S'^p{a., a'; t) 



yp{t)a]am, 

^aip{t)a\^, (38) 



with which we arrive at an operator form of the master equation shown below, 
p{t)^~'j^[H'{t),p{t)] 



r'(^) XI [2afcP(i)4' - alak'p{t) - p{t)alak-] 

1 

r(i) X[2afcp(t)4 - a\akp{t) - p(t)4afc], (39) 



fc=i 

where 



is the modified Hamiltonian of the two optical modes. From Eq. ([39]), we can see that 
besides the spontaneous decay of the individual mode, the environment, even only the vacuum 
fluctuation is considered, will result in a coherent interaction and a correlated spontaneous 
decay between the two modes. More importantly, our derivation of the master equation is 
fully non-perturbative, which goes beyond the Born-Markov approximation and contains all 
the back-actions between environment and the optical modes. The non-Markovian character 
resides in the time-dependent coefficients of the exact master equation. These formulae have 
been used to study the non-Markovian entanglement dynamics of two squeezed states j32j . 

The time-dependent coefficients in the exact master equation, determined by Eq. (j35p . 
crucially depend on the so-called spectral density, which characterizes the coupling strength 
of the environment to the system with respect to the frequencies of the environment. It is 
defined as J(lo) = \g-k\ ^(lo — uji). In the continuum limit the spectral density may have 
the form 

J{uj)^7jLj( — ) e--, (41) 

where ll>c is an exponential cutoff frequency, and 77 is a dimensionless coupling constant. The 
environment is classified as Ohmic if n = 1, sub-Ohmic if < n < 1, and super-Ohmic for 
> 1 [IS]- Different spectral densities manifest different non-Markovian dynamics. 
It is worth mentioning that an exact master equation has also been obtained very recently 
for the system of two harmonic oscillators bilinearly coupling with a thermal environment [31j , 
where the master equation is derived in the Wigner representation rather than the operator 
form of Eq. (|39p . Also the bilinear coupling in 31j is defined in terms of the coordinate 
variables of harmonic oscillators which is different from the coupling Hamiltonian we used 
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in Eq. ([4]). In terms of quantum optics language, the coupling between the system and the 
environment used in |31j involves simultaneously photon-photon scattering process and two- 
photon creation and annihilation process with the same coupling strength. Note that photon- 
photon scatterings are linear optical processes while two-photon creation and annihilation 
processes are non-linear optical processes, they cannot have the same coupling strength in 
quantum optics. Therefore, the model used in [31j might describe a physical system c^uite 
different from the optical system we considered in the present work. 

4.3 The Markovian approximation 

It is interesting to see that one can reproduce the conventional Markov solution from our exact 
non-Markovian master equation under certain approximation. By redefining the dynamical 
variables of the system as Q;/(r) = xi{T)e^^'^°'^ , and aJ(T) — X;(r)e*""'^, we can recast Eq. (|28l - 
[29| into 

XI+IKXV+ / dujJioj) / dT'e*('^«-")(^-^')[x, (r') 
Jq Jq 

+ Xxv (r')] = 0, (42) 

pOO pT 

i'i-iKx'i,+ / dujJ{uj) / dT'e-'^'^'>-'^'>^'-^">[x'i{T') 
Jo Jq 

+ \x[, (r')] = 0. (43) 

Then, we take the so-called Markov approximation, 

x{t')^x{t), x' {t') ^ x' {t) , (44) 

namely, approximately taking the dynamical variables to the ones that depend only on the 
present time so that any memory regarding the earlier time is ignored [52j . 

The Markov approximation is mainly based on the physical assumption that the correlation 
time of environment is very small compared with the typical time scale of system evolution. 
Also under this assumption we can extend the upper limit of the r' integration in Eqs. (|42j|43p 
to infinity and use the equality 



r — >oo 



lim / dr'e±*('^«~'^)(^-^') = 

T:5{uj-uo)Ti!^(—^), (45) 



where !^ and the delta-function denote the Cauchy principal value and the singularity, re- 
spectively. The integro-differential equations in (|42j|43p are thus reduced to a couple of linear 
ordinary differential equations. The solutions of xi and ij, as well as a/ and aj can then be 
easily obtained, which result in 



g — i(cJo— Ak)t _ g[-i(cJo+AK) — 2(7rJ((i)o)-i'5ii')]T 



2A 



(46) 
(47) 
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where 5uj = 3^ f^- 



3 u)—uq 



. Using the solutions (|46ll47p . one can verify from Eqs. p4)) that, 



T{t) -Ar'(t) - 7rJ(t^o), 

Q.{t) —uiQ ~ 5bj, fl'{t) — K — XSuj, 



(48) 



which is exactly the coefficients in the Markov master equation of the optical system f2? . This 
result can also be obtained easier by directly applying the Markov approximation Eqs. ()44ll45p 

to Eqs. jalETl). 

As shown above, all the coefficients in the master equation have become time-independent, 
and the non-Markovian master equation (|39p is reduced to the Markov master equation under 
the Markov approximation. This Markov approximation is valid to all kinds of spectral 
densities, including Ohmic, super-Ohmic and sub-Ohmic cases, while a different spectral 
density does produce the frequency shift, Suj — ^ , and decay rate, T = ttJ^loo), 

differently. As a result, we conclude that our exact non-Markovian master equation can not 
only explore more complicated situation where Markov approximation is unreachable, but 
also examine different spectral densities between the system and the environment even in the 
Markovian limit. This actually provides a simple way to reveal the underlying mechanism of 
quantum decoherence. 

5 Decoherence dynamics of entangled coherent states 

There are two different types of continuous variable entangled states. One is the entangled 
squeezed states, and the other is the entangled coherent states[5]. We have used the exact 
non-Markovian master equation derived in this paper to study the non-Markovian entangle- 
ment dynamics of two squeezed states in a recent published paper [32] • In this section, we 
will analyze the decoherence properties of the entangled coherent states. The decoherence 
dynamics of the two coherent modes is also fully described by the master equation ([39]) with 
the non-Markovian character residing in its time-dependent coefficients. The time-dependent 
coefficients in the master equation are determined by F± (t) as the solution of Eq. ([35]) for 
a specific environmental spectral density. In the present work, we will consider the Ohmic 
spectral density, i.e., n = 1 in Eq. (|4ip . which is often the case for optical communication. 

In Figs. [T] and [21 we plot the numerical results of the frequency shift Su!{t) and decay 
rate T(t) of the individual optical field as well as their corresponding Markovian values. It 
shows that the non-Markovian dissipation-noise dynamics is characterized by two time scales: 
Ti = 1/uJc (the shortest time scale of the environment) and T2 — I/ojq (the time scale of the 
optical modes). When t < ti, both coefficients, Suj{t) and T{t), grow very quickly, while after 
Ti, 5uj{t) and r{t) approach to the corresponding Markov values, given by Eq. ([48]) . gradually 
when the time approaches to the time scale T2. It clearly evidences that the non-Markovian 
effect has a huge deviation from the Markov effect within the time scale T2. This deviation 
will influence the dynamics later on significantly as a historical memory effect. The time 
dependent coefficients in the exact master equation (|39l) contain all the back-action effects 
between the system and the environment. The non-Markovian decoherence dynamics of the 
quantum optical field system thus becomes transparent due to the sensitive time dependence 
of these coefficients within the time scale T2 . 

In the following, we will investigate the decoherence dynamics of the entangled coherent 
states under the influence of the vacuum fluctuation. The entangled coherent states are 
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defined as 



\tp±) = l^ da, -a) ± I - a, a)), 



|<^±)=^=(|a,a)±|-a,-a)), (49) 



which were studied as quasi- Bell states |il[S3], where N± = 2(e2|"l ± g-^l"! ) are the nor- 
malization constants. Many schemes to generate such states have been proposed in optical 
systems and also in other systems |71 |S1 IHl HI] • It has also proposed to use these entangled 
coherent states for teleporting the superposed coherent states [TTl |2S] . 
The time evolutions of these entangled coherent states are given by 

^[e2(H^H''+(*)l^)(|a+(t),-a+(0)(a+(i),-a+(t)| + |-a+(t),a+(t))(-a+(t),a+(i)|) 
±g-2(|o|^-|a+(t)P) -a+{t)){-a+{t), -a+{t)\ + | - a+{t) , a+{t)) {a+{t) , -a+(t)|)" , 

±e-2(l"l'-l--(t)l')(|a_(t),a_(t))(_a_(t),_a_(t)| + | _a_(t),_a_(t)}(a_(i),a_(<)|)j50) 

respectively, where a±{t) = aF±{t). Eqs. ([50|) can be obtained directly from the exact solution 
of the reduced density matrix, Eq. ([TS]) plus Eq. ([?T|) by integrating over the initial variables. 
From Eq. ([55]) one can verify that for A = 1, the entangled coherent states \tp±) remain in 
pure states (decoherence free states) because a+{t) — aF+{t) — ae"*'^'^""'*^*, namely the time 
evolution of \ijj±) is independent of the decay rate T{t) and the shift frequency Suj(t) which are 
determined by F_ {t) when A = 1 [see Eqs. ([M]) and ([57)1 ] and only affect on the time evolution 
of the other two entangled coherent state \4>±). Similarly, when A = — 1, the entangled 
coherent states \4>±) becomes decoherence free states since a-{t) — aF^{t) = ae~*'-'^°+'*-'*, 
while the decay rate T{t) and the shift frequency 5u}{t) are determined by F+{t) which only 
affects the states \ip±)- Since A = 1 (or —1) corresponds to the case of the two optical modes 
coupling to the environment in phase (or out of phase), the above result indicates that two 
of the four entangled coherent states in Eq. (|49|) become decoherence-free states [15] if the 
two optical modes couple to the environment in phase (a constructive interference coupling) 
or out of phase (a destructive interference coupling). 

The reason that the V'-type and (p-type entangled coherent states in Eq. (j49)) have different 
decoherence behaviors comes from different symmetric properties of these entangled coherent 
states. The T/^-type and 0-type coherent states correspond to the center-of-mass and relative 
motions of two-field coherent states, respectively. This property becomes clear by defining the 
center-of-mass and relative motional variables of the two subsystems as = {a\ + a\ ) and 

= {a\ —02)- one can find, \^l'±) consist of only the relative motion, while \(j)±) lie only on 
the center-of-mass motion. When the two optical modes couple to the environment in phase, 
namely, gik = g2k = the interaction between the optical modes and the environment only 
affects the center-of-mass motion so that the entangled coherent states of the relative motion, 
\ip±), become decoherence-free states. On the other hand, if the two optical modes couple to 
the environment out of phase, i.e. gik = —g2k = 9k, the interaction between them only affects 
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the relative motion but leaves the entangled coherent states of the center-of-mass motion, 
free from decoherence. This is indeed a consequence of the sufficient condition for the 
decoherence-free space protected by symmetry [56]. 
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Fig. 1. Comparison of the decay rate V[t) [= Ar'(t)] between the non-Markovian (soUd hne) and 
Markovian (dashed hne) results. The parameters k/loq = 0.5, lJc/u>o = 30.0, and rj = 0.005 used 
in the numerical calculation. 
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Fig. 2. Comparison of the frequency shift (5a;(t) between the non-Markovian (solid line) and 
Markovian (dashed line) results. The parameters used in the numerical calculation arc the same 
as that in Fig. [l] 

We shall quantify the entanglement degree of the entangled coherent states by the familiar 
concept of concurrence usually used in a discrete basis [37]. To do so, we may rewrite Eq. ([50)1 
in terms of the orthogonal basis [53] , 



|0) =e 
II 



l«±W), 



h±(t)r 



(51) 
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with p±{t) — e^^l"*'-*-'! . Above change from the coherent state basis to the |0) and |1) 
basis is equivalent to a local unitary transformation of the states, which does not modify 
the entanglement degree in the original states. In this discrete basis, the concurrence can be 
calculated as usual. It is not difhcult to find that the concurrence (0) = (0) = 1, which 
is maximally entangled irrespective of the amplitude a, while C0^(O) = C^^(O) = tanh2 |a|^, 
which imply that the cf)^ and -0+ states are not initially maximally entangled. One can also 
show that when the two optical modes couple with the environment in phase, i.e. A = 1, the 
concurrence (t) = 1, and C^^.(t) = tanh2 |q;| , namely, the entanglement of \'ip±) remain 
unchanged during the time evolution. While \(l)±) are sensitive to decoherence. In contrast, if 
the two optical modes interact with the environment out of phase, i.e. A = — 1, (t) ~ 1, 
and Cip^{t) = tanh2 |a|^, while C^_(t) and C^^(i) will decay (disentanglement) due to the 
decoherence. 

In Fig. 3, we show the concurrence evolution in time for the entangled coherent states 
and IV'i). With the in-phase coupling between the optical modes and the environment 
(A = 1), our numerical results verify that the entanglement degrees of \(j)±) (given by the 
solid and dot-dashed lines in Fig. [3]) suffer from a fast decay during the time evolution while 
the entanglement degrees of jV'i) remain unchanged (the dot-dot-dashed and dot-dot-dot- 
dashed lines in Fig. [3|). To compare the non-Markovian entanglement dynamics with the 
Markovian dynamics, we also plot the concurrence evolution for \(l)±) under the Markovian 
approximation, denoted by the dashed and dotted lines, respectively, in Fig. [3l As one can 
see, the non-Markovian effect accelerates the disentanglement. This is mainly a contribution 
of the short time peak in the decay rate r{t) as a memory effect. It is also worth noting that 
no entanglement oscillator is observed in above solution even there has coherent coupling fl'(t) 
presented. This is because the two-optical-field coupling n'{t) contributes only a global phase 
to the entangled coherent states during the time evolution, which has no influence on the 
entanglement degree of the states. While as expected, |'0±) are decoherence-free irrespective 
of Markovian or non-Markovian dynamics being considered. 

For the case of out-of-phase coupling between the optical modes and the environment, 
namely A = — 1, the roles of the decoherence effect on the 0-type and the f/'-type entangled 
coherent states are exchanged. The (j)-type entangled coherent states remain unchanged, while 
the tp -type entangled coherent states are disentangled by decoherence. The numerical results 
of the entanglement evolution for \'ip±) and \(j)±) are given by the same curves in Fig. [3] with 
the exchange between ji/ii) and \(j)±) states as in the in-phase coupling case. 

From the above analysis, one can find that when the two identical optical modes couple to a 
common environment in an arbitrary phase difference, no decoherence-free entangled coherent 
state can exist among the four entangled coherent states. All the four entangled states could 
be disentangled by the vacuum fluctuation in time. The non-Markovian dynamics will speed 
up the disentanglement process with respect to the Markovian approximation. However, the 
parameter A models the phase difference between the two optical modes coupling with the 
environment. Physically, it is always possible to adjust the two optical modes such that the 
couplings of the two optical fields with the environment are either in phase (A = — 1) or out 
of phase (A = —1). Then two decoherence free states among the four entangled coherent 
states can always be constructed in principle. It is certainly interested in seeing experimental 
evidences on the preservation of two decoherence free entangled coherent states as well as the 
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Fig. 3. Time evolution of the concurrences for different initial states. The solid and dot-dashed 
lines show the non-Markovian time evolution of the concurrences for 10—) and \4>+) with A = 1 (or 
and with A = —1), respectively. Their corresponding Markovian behaviors are shown 

as the dashed and dotted lines, respectively. The dot-dot-dashed and dot-dot-dot-dashed lines are 
the concurrences for \il>—) and with A = 1 (or \4>~) and \4>+) with A = —1), respectively, 

which remain unchanged during the time evolution. The initial coherent state parameter a = 0.8, 
and other parameters are the same as in Fig. [Tl 

non-Markovian disentanglement enhancement to the other two entangled coherent states. 
6 Summary and Discussions 

In summary, we have studied the detrimental effects of environment on the entangled coherent 
states. We microscopically modeled the decoherence dynamics of entangled coherent states 
under the influence of vacuum fluctuation. An exact master equation with time-dependent 
coefScients reflecting the full memory effect of the reduced system has been derived by using 
the Feynman- Vernon influence functional theory in the coherent-state path-integral represen- 
tation, which enables us to treat both of the back-actions from the environment to the system 
and from the system to the environment self-consistcntly. In addition, we have also explicitly 
deduced to the well-known Markovian dynamics for the optical modes from our exact non- 
Markovian master equation in the Markovian approximation. The analytical analysis of the 
difference between the non-Markovian dynamics and its Markovian approximation presented 
in this paper may provide a quantitative way to experimentally explore the non-Markovian 
effect as well as the spectral densities between the system and the environment even in the 
Markovian limit. 

We then investigated the non-Markovian dynamics of the entangled coherent states, one 
of two typical continuous variable entanglement states often used in quantum information 
processing. The other type of continuous variable squeezed states has already be studied 
by two of us based the same master equation derived here [52] ■ Our first-principle analysis 
shows that the non-Markovian effect accelerates the disentanglement compared with the re- 
sults based on Markovian approximation. It is the short time peak of the time dependent 
coefficients in the master equation, which is incorporated with the system's dynamics as a 
historical memory effect, that contributed to this acceleration. Although the Born-Markovian 
approximation has been widely employed in the field of quantum optics, we argue that our 
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investigation might be helpful for understanding decoherence in nanoscale cavity devices and 
ultrafast optical processes. For example, we have noticed the rapid development of optical 
cavity technology, which has been employed to confine a single atom [57] or a single quan- 
tum dot [Ml [59] in strong coupling regime, the prerequisite of quantum network. The strong 
interaction occuring in nanometer size subject to vacuum fluctuation suffers from some un- 
predictable incoherence errors [57], which should involve the non-Markovian effects. In this 
sense, our study of non-Markovian dynamics, although with a simplified model, paves a way 
toward clarification of the mechanism regarding those incoherence sources. 

We have also shown how the decoherence behaviors of the different entangled coherent 
states depend on the symmetrical properties of these entangled coherent states as well as 
the interference properties of couplings between the two optical modes with the vacuum elec- 
tromagnetic environment. Since the exact non-Markovian master equation has been derived 
non-perturbatively and exactly, decoherence dynamics subject to different spectral densities 
of environment would be naturally available by our treatment. In fact, the non-Markovian 
master equation (i.e., Eq. (|39p ) derived in this paper has been used in treating the decoherence 
dynamics of entanged squeezed states with sub-Ohmic and super-Ohmic spectral densities of 
the environment [32| . More complicated cases, e.g., the environment at finite temperature, 
would be hopefully figured out by the similar way to the derivation of Eq. ([3T|) . As a final 
remark, we would like to mention a very recent experiment for distinguishing different co- 
herent states |60j , which shows a potential of using entangled coherent states for quantum 
communication. The entangled coherent states have significantly different properties from 
the entangled squeezed states and have been proposed as another type of continuous variable 
quantum channels. As robustness of the quantum channel is essential in view of decoherence, 
we expect that our consideration of decoherence dynamics of entangled coherent states would 
be useful for understanding quantum communication experiments with realistic ultrafast op- 
tical processes in nanocavities and nanophotonic systems. 
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